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Abstract 

We report observations of the gyrotropic change in magnetic susceptibility of the 
Fe 3+ electron paramagnetic resonance at 12.037GHz (between spin states [1/2 > 
and 1 3/2 >) in sapphire with respect to applied magnetic field. Measurements 
were made by observing the response of the high-Q Whispering Gallery doublet 
(WGH±nfi,o) m a Hemex sapphire resonator cooled to 5 K. The doublets initially 
existed as standing waves at zero field and were transformed to traveling waves due 
to the gyrotropic response. 



*Electronic address: kbenmess@cyll ene.uwa.edu.au 
'''Electronic address: mike@physics.uwa.edu.au 



I. INTRODUCTION 



Residual paramagnetic impurities in ultra-pure Hemex sapphire crystals greatly influence 
the electromagnetic properties of the resonators cooled to near liquid helium temperature. 
For example, the electron paramagnetic resonance (EPR) of impurities such as Cr 3+ , Ti 3+ 
and Mo 3+ , influence the temperature dependence and allow frequency-temperat ure turnover 
(annulment) near l iquid helium temperature , which can vary betwe en 2 to 10 K (j Jones et al. 



1988 



Mann et al. 



1992 



Luiten et al. 



1996 



Kovacich et al 



, 119971 ) . This phenomenon is vi- 



tal to produce u l tra-stable frequencies for cryogenic sapph i re osc illators fjGiles et al 



Bourgeois et al. 



2004; 



Hartnett et al. 



2006 



Locke et al. 



1990 



20081 ) an d the study o f EPR 



in low-loss c r ystals has become an important topic of investigation (jHartnett et al. 



1999 



Luiten et al. 



2001 



Hartnett et al. 



2007 



20011 ) . More recently, it was shown that the small 



amount of residual Fe 3+ ions (less than parts per million) are sufficient to create Maser 
oscilla tions due to the zero-field 3- level system between j l/2 >, 1 3/2 > and 1 5/2 > spin 



statesl jBourgeois et al 



2005 



2006 



Benmessai et al. 



20071 ). This has allowed a new way to 



operate a cryogenic sapphire oscillator, where population inversion is created by pumping 
the spin states in the 1 1/2 > transition to the 1 5/2 > transitions with a 31.3 GHz pump, and 
has the po tential to operate with a frequency stability governed by the Schawlow-Townes 
noise limit (jBenmessai et a/1120071 ) (sub 10 -16 frequency stability). 

In this work we show how the application of an axial magnetic field on a cryogenic sapphire 
resonator, with a mode tuned on the Fe 3+ EPR, adds a gyrotropic component of magnetic 
susceptibility. This is achieved due to the very high-Q whispering gallery (WG) modes (of 
order 10 9 ) that are excited, which enables the discrimination of the two degenerate standing 
wave modes known as a doublet. Precision measurement of the frequency and Q-factor of the 
doublet reveals an asymmetric response, which demonstrates that the susceptibility added 
by the magnetic field is predominately of gyrotropic origin and that the modes transform to 
traveling waves. Traveling waves are known to enhance the interaction between the pump 
radiation and the maser medium due to the elimination of standing wave nodes, and should 
be a way of generating greater output power in the respective solid-state maser. 
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II. WHISPERING GALLERY MODE FIELDS 



To determine the electromagnetic fields of the mode under investigation we use a sep- 
aration of variables technique, which has proved to be reliable for the approximate calcu- 



lation of e 



ectromagnet i c fie 



resonators ( Tobar et al. 



1991 



ds and properties of high-Q WG modes in low loss sapphire 



Wolf et all |200J). To calculate the gyrotropic properties we 
apply perturbation analysis on the non-gyrotropic fields. The technique assumes the res- 
onator is a perfect cylinder of uniaxial anisotropy with the c-axis aligned along the cylinder 
z-axis as shown in Fig. [TJ 

For electromagnetic modes exisiting in non-gyrotropic media with one or more azimuthal 
variations (m > 0), two degenerate mode solutions exist (doublet) as orthogonal displaced 
standing wave solutions proportional to either cos(m<j)) or sm(m0). The degeneracy of the 
modes is lifted by a very small amount (only observable in high-Q systems) due to dielec- 
tric backscatter or other perturbations, which couple the counter clockwise ( CCW, e jm ' 



or clockwise (CW, e" 



-jm4 



1995 



Kippenberg et al. 



parts between the two double t standing wave modes ljWeiss et al. 



2002 



Gorodetsky et al 



2000). In gyrotropic media, the doublets 



have their degeneracy further modified in a non-reciprocal fashion, and the modes must 
transform to the travelling wave basis (rather than standing wave) with o pposite ellipti 



calpo 



1996.? 



arization due to the non-reciprocity imposed by the gyrotropic effect (IKrupka et al. 



In this work we observe gyrotropic effects on a quasi- Transverse Magnetic (TM) WG 
doublet pair, which is represented by the notation, WGH ±m ^^. Here n the number of 
radial zero crossings and p the number of axial zero crossings and the '±' symbol represents 
the sense of the elliptical polarization of the doublets. For the fundamental WGH± mi o,o 
mode family the E z electric field is the dominant component and is symmetric along the 
z-axis (around z = 0) . Given these conditions one can show that the solutions to Maxwell's 
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(1) 



equations yield (the time dependence e?" 1 is assumed): 

E zl± = A l± e ±jm *cos(f3z)J m (k E r) 
E z2± = A 2± e ±jm *cos(Pz)K m (k out r) 
E z3± = A^e^e-^JUkEr) 
H zl± = B l± e^ m4> sin^z)J m {k H r) 
H z2 ± = B 2± e ± ' m<t 'stn((3z)K m (koutr) 
H z3± = B^e^e-^JUknr) 

To calculate the other components we apply the following Maxwell's relationships to the 
z-components given in Eq. (CQ) in the regions 1, 2 and 3 (derived from Maxwell's equations 
in cylindrical co-ordinates). 



k 2 + 



k 2 + 



dz 2 

d 2 



d 2 H : 



H r ± = juJ6 ± e -^^ + V?* (2a) 
r ocp ozor 

dH z± 1 d 2 E z± 

E<t>± = J^O^— + ~ T^Zal (2b) 

(2c) 

//- + ^ ) E r± = -jcufio^^ + (2d) 
oz z J or r ozor 

Here, en = 11.349 and e± = 9.272, are the uniaxial permittivity components of sapphire 

cooled close to liquid helium temperature (~5 K) parallel and perpendicular to the c-axis 

respectively. 

The mode pair that we are interested in is the WGH±i 7i0 ,o inside a sapphire crystal of 50 
mm diameter and 30 mm height at room temperature. Contraction from room temperature 



dz 2 



Or 



r dzdcj) 



Q 2 



d 2 



-jue±e — 
r 



1 dE 7+ d 2 H z± 



+ 



dcj) dzdcj) 



dH z± ld 2 E z ± 



to 5 K is taken into account in the modelling by the fac t ors of 0.99927 along t 



0.99939 perpendicular to the z-axis (lTobar et al. 



1997 



Krupka et al 



1996 



re z-axis and 
Matching 



20041 1 . we obtain 



boundary conditions with the same approximations as given in (1 Wolf et al. 
the following values of the parameters as shown in Table [H note that initially no gyrotropic 
effects are assumed (the effects are small and perturbation analysis suffices). This means 
the initial calculation of the + and - polarized modes are degenerate. 

The parameters in Table [I] are independent of the elliptical polarization, and electromag- 
netic filling factors may be calculated from Eq. (JHJ), which also turn out to be independent 
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of elliptical polarization, and are shown in Table HT1 



Pr, 



jjj v ^H-Hdu V - JJJ v e(u)E.Edu 

Here % defines the component of the field and corresponding component of filling factor 
within the corresponding region of volume VJ. The normalization factor calculates the total 
field energy across all regions in the total volume of the resonator, V. 

Note that the largest components of the field within the sapphire are E z and H r , this 
means that the dominant propagation is around the azimuth (Poynting Theorem) as ex- 
pected for a WG mode. Also, the three most dominant components (H r , and E z ), are 
that of a TM mode with the component combining with the H r to define the sense of the 
elliptical polarization. The amount of field in the vacuum surrounding the dielectric (regions 
2 and 3) may be calculate by, 1 — p ez — p er — p e ^ = 0.00679954, and 1 — p mz — p mr — p m ^ = 
0.0828332 respectively. 

To determine the sense of elliptical polarization we take a closer look at the solution of 
the dominant components for the + and - mode in region 1 (the sapphire dielectric), which 
are calculated by substituting Eq. (CQ) into (J2J): 



(3) 



E zl+ = e^cos(Pz) \-E zlo J m (k E r)] 
H rl+ = e jm *cos((3z) 



dJ m (k H r) J m {k E r) 

Hrla o r H r lb 



H^ 1+ = je^costfz 



H, 



01a ' 



dr 

dJ m {k H r) 
dr 



016" 



r 

J m {k E r) 



E z i- 
H r i- 

H(l>i- 



-jm<t> 



cos((3z) [E zlo J m (k H r) 



(4) 



-jm<t> 



cos(j3z) 



dJ m (k H r) Jm{k E r) 

tlrla o r H r lb 



-jrrvjs 



cos((3z) 



dr 

dJ m (k H r) 
dr 



016" 



r 

J m {k E r) 



Here, the solution requires H rla = 1.09 x 10 E zlo , if, 



ib 



r 

1-4, 



1.82 x 10- 4 ^ lo , H&a = 1.86 x 
10~ 6 -E z io, H^ii, = 1.07 x l0~ 5 E z i o . Thus, by inspection the polarization of the magnetic field 
for the CCW propagating mode is proportional to e 7 "^ and has CCW elliptical polarization 
(+ mode), while the magnetic field for the CW propagating mode proportional to e~^ mr ^ has 
CW elliptical polarization (- mode). 
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III. GYROTROPIC REPRESENTATION 



In this work we show that gyrotropic effects occur on the WGH 17,0,0 doublet pair when 
tuned within the Fe 3+ electron paramagnetic resonance in the presence of a DC mag- 
netic fiel d. To calculate the gy rotropic nature we use the same notation as presented by 



Gurevicht lGurevich et al. 



19601 ) . We assume the case of a weak magnetic response, which is 



true for our situation as we use the purest HEMEX sapphi re available, which has no more 



than t he order of parts per million of residual im purity ions (IJones et al. 



1988 



Mann et al. 



1992 



Hartnett et al. 



2007 



Tobar et al 



20031 ). In this case the magnetic susceptibility 



added can be considered independent of the microwave field, i.e. only dependent on the 
applied DC field. The permeability added by the electron paramagnetic resonance response 
to the magnetic field in the axial direction, is a tensor of the form: 



fi ± (B) jyL a {B) 

-3lia{B) (JL ± (B) 

1 



(5) 



Thus, the change in permeability (or change in susceptibility) created by the response to 
the external field is given by the following: 



AV( J B) = V( J B)- y u / 



A/i ± (5) jfi a {B) 
-jfx a (B) Afi ± (B) 




(6) 



Here, I is the unitary matrix and in general the permeability components are complex. 



IV. APPLICATION OF PERTURBATION THEORY 

The typical value of susceptibility added by residual impurities is of the order 10~ 8 , thus 
we apply perturbation theory based on the fields calculated using separation of variables 
in Sec IHII Following the method in Gurevich, the change in complex frequency maybe 
calculated from the change in complex permeability from, 

A/ ± = JJJ Vl H±A*J? H±dis 

/± fffy e{v)E ± ■ E ± dv + jjjy f, H ± ■ H ± dv [ } 
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Here, fo is the unperturbed frequency (no gyrotropic effect) so that A/ ± = f ± — /o, where 
Af ± /f ± is in general complex. In this work the complex fractional frequency is defined by 
A/± / f± = Af Re± I f ± +jAf Im± I f ± , where the real part describes the frequency shift and the 
imaginary part the unloaded Q shift, where Qg^ ~ 2Afj m± /f ± . The total mode unloaded 
Q values (Qo±) is related to the zero field unloaded Q value (Qzf±) an d t ne shift due to the 
magnetic field by Qo^=Qzf ± +Qb1- For the above model the complex permeability added 
by the magnetic field is of the form A[i± = A/i' ± — jA/i" and \i a = f/ a —j/jf^. Substituting the 
complex permeability into Eq. ([6]), and then ([6]) into ([7]) along with the unperturbed fields, 
the following sensitivities with respect to the permeability are calculated for the WGH-1-17 0,0 
doublet. 

A/ fie+ 



— - -0.4571 76 (A// ± - 0.10485*0 = -0.457176A/4 //+ 

A ^ Re ~ = -0.4571 76 (Afi' ± + 0.10485*0 = -0.457176A// e// _ 

Aflmj- 



0.457176(A//[ - 0.10485/O = 0.457176 A/i" //+ 



(8) 



/+ ^ ^ ™ r~e„ + 2Qb _ 

A fi™- _ o.457176(A//[ + 0.10485^') = 0.4571 76 A^ ff _ - 1 



The non-reciprocity of the two modes is highlighted by Eq. (jSJ), which shows the - and + 
mode dependent on different effective permeability A// e // ± = A/0^ ± — j A/x^ ± (both real 
and imaginary). 

The magnetization model of Gurevich predicts n' ± and \J' a to be of the same sign. Thus, 
given that the both modes tune up in frequency with magnetic field the '-' or CW mode must 
be the mode with the greatest tuning coefficient. According to the model the mode should 
also be accompanied with a decreasing Q-factor, while the '+' or CCW mode increases its 
Q-factor. This is indeed what is observed, in the next section two resonators are analysed 
and the values of A/i' ± , /jf a , A{i' ± and fi' a calculated as a function of magnetic field. 



V. EXPERIMENTAL RESULTS 



Hemex sapphire crystals are known to have Fe 2+ impurities in the crys tal. To convert 



these ions to Fe 3+ the crystals need to be annealed in an oxygen environment (IBenabid et al. 



20001 ) . The resulting ESR exhibits a zero fiel d ESR between the |l/2 > and 



at 12.037 GHz with a 27 MHz bandwidth tjBourgeois et al 



2005 



2006 



3/2 > spin states 



Benmessai et al. 



20071 ). The crystal resonators are nominally 50 mm diameter and 30 mm in height, for 
such a resonator the WGH±i 7i0) o doublet pair is tuned inside the ESR bandwidth, with a 
predicted frequency from the separation of variables technique of 12.03 GHz. In this work we 
report on the properties of two crystals, which were annealed. Measurements were taken on 
a crystal, which was not annealed, however no measurable effects with respect to magnetic 
field were observed. The experimental set up of the cryogenically cooled system with a 
variable external magnetic field is shown in Fig. [2j 

The Q-factors and frequencies of the two WGH ±17 0j0 doublet modes were recorded as 
a function of magnetic field for two separate sapphire resonators using a network analyser 
in transmission. A typical recording of the transmission coefficient with and without the 
magnetic field on is shown in Fig. [3J 

With zero field applied to the crystals the following parameters in Tab. 3 and 4 were mea- 
sured. The doublet mode s were already separated in frequency at zero-field due to other non 



gvrot ropic perturbations llKippenberg et all 120021 ; iGorodetsky et al\ . l200ol : iBourgeois et al. 



20051 ). The modes are labelled as upper (subscript up) or lower (subscript low) for the higher 
and lower frequency tuned modes. The frequency shift and unloaded Q-factor of the upper 
and lower modes were measured as a function of magnetic field and are shown in Fig. H] and 
[5] for the two resonators. 

From the properties present in Fig. H] and it is straightforward to determine the 
polarizations of the modes. The '-' CW polarized mode has a larger effective permeability 
with respect to magnetic field and should thus tune at a greater rate. The same mode 
should also increase its effective loss tangent at the same time. Thus, for sapphire No. 1 
the '-' mode is identified as the upper mode and the '+' CCW mode as the lower. It is the 
reverse situation for sapphire No. 2 as the lower mode is the '-' mode and the upper is the 
'+'. After this identification, equation [9]) may be used to solve for the real and imaginary 



permeabilities. The solutions lead to the following equations: 

' S7 \7^— + 7^—) 

(9) 



A//' : -.1.09307 ( + ) A// j — 1.09367 ( -— ^ h j^— 



A/ Rc - 

f- 






A/ Rc - 


A/ Re+ \ 




/- 


/+ J 





fi' a = -10.4305 -^P ^ = 10.4305 



1 1 



,Qb- Qb- 

Thus, from Eq. and the data of Fig. H] and 0, the change in complex permeability with 
respect to magnetic field may be calculated and is plotted in Fig. [6] and [71 

It is clear that the gyrotropic susceptibility added by the magnetic field is larger than 
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the non-gyrotropic component for both real and imaginary terms. 



VI. TRAVELLING WAVES 



Coupled mode theory is used to show that standing waves are converted to travel- 
ling waves. Previous analysis has shown that in high-Q systems perturbations, including 



backscatter due to the dielectric and other p erturbations such as probes, and wil 



the WG doublet to have its degeneracy lifted (j Weiss et al 



Mazzei et al. 



1995 : 



Gorodetsky et al. 



cause 



2000 : 



20071 ). Co upled modes in tuneable h igh-Q sapphire resonators have been stud- 



ied in detail previously ljTobar et al. 



1991 



19931 1 . The lifting of the degeneracy is caused 



by reactive coupling through the electromagnetic fields (strong coupling), which creates a 
band gap between the modes. The frequency of the upper and lower modes with respect to 
a linear tuning parameter, x, are in general given by Eq. fflOl) . 

fu P /ioM = 1 ± ^ i 

(10) 



fo 



/+(*) 



1 + kx 



/-(*) 



f A 2 /4 
1 — kx 



fo fo 

Here the frequencies of the travelling waves /+ and /_ depend linearly on the tuning pa- 
rameter x, with a tuning coefficient of ±k respectively, an idealized diagram with respect to 
coupled mode theory with k = 1/2, is shown in Fig. [HJ When the two modes have the same 
frequency, i.e. x = the mode separation due to the strong coupling is f up / fo~ flow/ fo = A, 
which defines the strength of the coupling and the band gap width. 

The data represented in Fig. H] and [5] shows the mode frequencies tuning as a func- 
tion of magnetic field. The magnetic field tunes the frequencies due to the tuning of the 
EPR, however the tuning is non-linear. Given that the non-gyrotropic response of the 



EPR to the magnetic field s 



r ould be in the linear regime, as the magnetic fields are rela- 



tively small (jKornienko et al. 



19611 ) we may define the average frequency shift of the two 



modes as an effective linear tuning parameter, such that x = B/\B\(f up + fiow)/2fo- Fig. 
M and [TU] show the differential shift between the two modes as a function of x (ignoring 
the initial mode splitting). The differential shift according to Eq. (fTOl) will be of the form 
(fup — fiow)/2fo = \/k 2 x 2 + A 2 /4, and may be fitted to the experimental data to calculate 
the coupling parameter, A, between the travelling waves. 

The ESR is inhomogenously broadened so the interaction will be with the spin packet of 



closest frequency in both sapphires. Sapphire 2 has a lower threshold of about 0.6 Gauss 
compared to sapphire 1 with 4 Gauss. The threshold value basically corresponds to the 
value of magnetic field that the standing waves are transformed to travelling waves and the 
point where the mode responds linearly to tuning parameter x with the value of k for the 
tuning coefficient. The threshold magnetic field value also correlates with the value of the 
coupling strength between the travelling waves with the higher value demanding a stronger 
coupling. Thus, the more tightly coupled the travelling waves the stronger the magnetic 
field threshold that is required to transform them into travelling waves. 

The doublet splitting at zero-field is of order 10~ 7 , which is more than two orders of 
magnitude than the coupling parameter measured in this work. However, this splitting is 
most likely due to the non-perfect cylindrical radius. Once the modes transform to standing 
waves they sample orthogonal space along the azimuthal dimension and will "see" a different 
effective radius. Only a difference of 2.5 /xm over the 25 mm radius is required to account for 
this. The hypothesis is supported by the fact that the tolerance of the radial circumference 
in the machining process given by the manufacture is ±50 /xm. 

VII. CONCLUSION 

The contribution to the magnetic susceptibility (or permeability) of the electron paramag- 
netic resonance of Fe 3+ ions in sapphire at 12.03 GHz (between spin states 1 1/2 > to |3/2 >) 
with respect to a DC axial magnetic field has been measured from precision measurements of 
Q-factor and frequency of a whispering gallery mode doublet pair at 4.2 K. The form of the 
susceptibility was shown to be predominantly gyrotropic causing non-reciprocal response 
of the measured resonant modes. We showed that the gyrotropic response converted the 
standing waves to travelling waves due to the non-reciprocal response. 
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FIG. 1 Solution regions for separation of variables in cylindrical coordinates. 
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Gold platted copper cavity 
Anti-radiation shield 

FIG. 2 The sapphire resonator was mounted in a gold platted copper cavity and measure in 
transmission with 10 dBm of incident power from a network analyser. The DC magnetic field in 
the axial direction was applied by a Helmholtz coil pair system as shown above. The resonator 
was placed in vacuum and fixed on a cryocooler (PT 405) cold finger and cooled close to 5 K using 
a temperature control system. 
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Frequency (GHz) 

FIG. 3 Frequency, Q factor and mode coupling were recorded using a network analyser referenced 
to a Hydrogen maser to avoid frequency drift of the network analyser. For resonator 1 the coils 
sensitivity is 4 gauss/ A and 1 gauss/A for resonator 2. 




B (Gauss) B (Gauss) 

FIG. 4 Left, Measured fractional frequency shift for the WGH-1-17^,0 lower and upper doublet modes 
as a function of applied axial magnetic field for sapphire No. 1 (Table 3). Right, corresponding 
measured unloaded Q-factor. 
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B (Gauss) B (Gauss) 



FIG. 5 Left, Measured fractional frequency shift for the WGH-ti7 o,q lower and upper doublet modes 
as a function of applied axial magnetic field for sapphire No. 2 (Table 4). Right, corresponding 
measured unloaded Q-factor. 




B (Gauss) B (Gauss) 

FIG. 6 Left, Real permeability components as a function of applied axial magnetic field for sapphire 
No. 1. Right, corresponding imaginary components. Note, the gyrotropic response is larger in 
both cases. 
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FIG. 7 Left, Real permeability components as a function of applied axial magnetic field for sapphire 
No. 2. Right, corresponding imaginary components. Note, the gyrotropic response is larger in 
both cases. 




-0.4 -0.2 0.2 0.4 

FIG. 8 Normal mode frequencies versus tuning parameter, x, with k = 1/2 and A = 0.1 as given 
by Eq. ([TO]) . The top solid curve is the upper tuned mode f up /fo, the bottom solid curve is the 
lower tuned mode fiow/fo- The non-coupled modes in the travelling wave basis / ± //o are given 
by the dashed lines. When the modes are detuned they behave as travelling waves, when they are 
tuned the modes interact and are mixed and are best described by standing waves (cos(m(f>) and 
sin(m(f))). The extent of the transition region depends on the strength of the coupling. 
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FIG. 9 Fractional differential mode frequency shift versus common mode frequency shift for sap- 
phire No. 1. The curve fit gives the coupling parameter (m2) of A = 7.8(1.9) x 10 -10 and coefficient 
(ml) k = 0.95(0.05). 
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FIG. 10 Fractional differential mode frequency shift versus common mode frequency shift for 
sapphire No. 2. The curve fit gives the coupling parameter (m2) of A = 1.5(2.4) x 1CT 11 and 
coefficient (ml) k = 0.40(0.01). 



Tables 

TABLE I Calculated parameters from matching boundary conditions for the WGH-1-17^,0 mode 
pair. 

Frequencyi kn k E (3 k out an «£ 

12.0305 GHz 717.707 802.952 103.866 760.709 841.609 j 229.754 

TABLE II Electromagnetic filling factors within sapphire dielectric (region 1) for the WGH±i7 ; o,o 
mode pair. 

Pezl Perl Pe</>1 Pmzl Pmrl Pm(j>l 

0.970906 0.00342936 0.0188647 0.0029268 0.835142 0.0790977 

TABLE III Zero field properties of the WGH±i7 ; o,o doublet in sapphire resonator number 1. 
Mode (Sapphire No.l) Frequency [GHz] Q 

lower 12.023984181 3.5 x 10 s 

upper 12.023994100 13.8 x 10 s 

TABLE IV Zero field properties of the WGH±i7 ; o,o doublet in sapphire resonator number 2. 
Mode (Sapphire No.2) Frequency [GHz] Q 

lower 12.041856348 5.9 x 10 s 

upper 12.041861773 7.3 x 10 8 
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